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Abstract. We consider the multi-particle lattice Anderson model with an i.i.d. 
random external potential and a short-range interaction. Using the multi-particle 
multiscale analysis (MPMSA) developed by Chulaevsky and Suhov [6], we prove 
spectral localization for such Hamiltonians at low energies under the assumption 
of log-Holder continuity of the marginal probability distribution of the random 
potential. Under a stronger assumption of Holder continuity, Anderson localization 
for such systems at low energies was established earlier by Aizenman and Warzel 
(2009) 1 with the help of the multi-particle Fractional-Moment Method. 



1. Introduction 

Multi-particle Anderson localization theory is a relatively recent direction in the 
spectral theory of random operators. The general structure of an ./V-particle Hamil- 
tonian of a lattice quantum system with interaction is as follows: 

N 

H m H = -A + V(. Xj ,u) + Ufa,. ..,x N ) (1.1) 

3=1 

where A is the nearest-neighbor lattice Laplacian in (Z d ) N ~ Z Nd , V : Z d x ft — > K is 
a random field on Z d relative to some probability space (SI, P), x%, . . . ,xn G Z d , and 
U: (Z d ) N -4 R is the potential energy of inter-particle interaction. 

The first mathematically rigorous results have been obtained by Chulaevsky and 
Suhov [51|B] with the help of the Multi-Scale Analysis (MSA) and by Aizenman and 
Warzel pQ who used the Fractional- Moment Method (FMM). In both cases, it was 
assumed that the random potential field is i.i.d., and the interaction U has finite 
range. In [6], the multi-particle spectral localization was proven for strongly disordered 
systems. The method of PQ[5] allows to establish dynamical localization in the following 
situations: 

• Strongly disordered systems — for all energies. 

• Weakly disordered systems — for "extreme" energies. 

• Multi-particle systems with weak interaction, under the assumption that the single- 
particle subsystems are localized in absence of interaction. 

Due to technical requirements of the FMM, it was assumed in [T] that the i.i.d. 
random field V has Holder-continuous common marginal distribution. 
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In the present paper, we show, using Multi-Particle Multi-Scale Analysis (MPMSA), 
that this assumption can be relaxed to log-Holder continuity of the marginal distribu- 
tion: see assumption (P2) below and assertion (B) of Theorem [TJ 

On the other hand, under a more restrictive condition — assumption (PI) given 
below, one can apply a simpler multi-particle scaling procedure based on a recent eigen- 
value concentration bound proven in [3] and establish exponential spectral localization 
along with more efficient intermediate bounds on the decay of eigenfunctions in finite 
cubes: see assertion (A) of Theorem Q] 

Our main results are Theorem [T] (multi-particle localization at low energies), Theo- 
rem [2] and Corollary [1] (stability of multi-particle localization for weak interaction). 

Section [5] is devoted to the description of the multi-particle model and to the main 
statements. Assertion (A) of Theorem[T]is proven under assumption (PI) in Sections|2] 
andUl The parts of Theorem [T] under assumption (P2) are proven in Section[S] Finally, 
Theorem [5] and Corollary Q] are proven in Section [Bl 

2. The model and the main results 

2.1. Basic notations. We are interested in the A^-particle hamiltonian (jl.ip under 
some assumptions on A, V and U. In what follows N > 2 is fixed. Since multiscale 
analysis applied to this hamiltonian requires also the consideration of hamiltonians for 
rt-particle subsystems for any 1 < n < N, we introduce notations for such n. 

2.1.1. The Laplacian in Z y . We generally equip 1 V with the max- norm | • | defined by 

\x\ = max-Oil, . . . , \x„\}, 

for x — (x\, . . . : Xv) 6 Z^, but, in the definition of the laplacian we will use another 
norm defined by 

Mi = |xi| H h \x v \. 

We denote by A the ^-dimensional lattice nearest-neighbor Laplacian acting on the 
Hilbert space £ 2 (Z U ) by 

(A¥)(z)= Yl - *(a0) = -2i/¥(a;) + J2 

|y-x|i=i |y-x|i=i 
Note that A is bounded and —A is nonnegative. 

2.1.2. The n-particle Hamiltonian. Let n > 1 and d > 1 be two integers. A configura- 
tion of n distinguishable quantum particles in the lattice Z d is represented by a lattice 
vector 

x= ( Xl ,...,x n ) € (Z d ) n =Z nd , 
where Xj — (x^ , . . . , Xj ) G Z d , j = 1, . . . , n. We will consider a random Hamiltonian 
H^™) of n particles in Z d of the form 

n 

HW(u) = -A + J]% 1 u) + U ) CU) 

J'=l 

acting in the Hilbert space = i 2 (Z dn ). It is easy to see that both norms are 

compatible with the identification Z dn ~ (Z d )™, i.e., 

x| = max \xj I, 

l<j<n 

|x|i = I a? 1 1 1 H h \x n \i. 
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Hence we have 

n 

A = V Ip (Z£i) <g> • • ■ <g> l P(zd) <g)A (j) ® Ip (Z£!) ® • • ■ <gi I £ 2 (Zt!) , 
3=1" " ' * < ' 

j — 1 times n— j times 

where A^ denotes the lattice nearest-neighbor Laplacian acting on the jth factor Z d 
of the decomposition (Z d )™ = Z rf x • • • x Z d , and Imgd) is the identity operator on 
J? 2 (Z d ). 

2.2. Basic assumptions. 

2.2.1. Assumption on A. 

(K) Kinetic energy. A is the 7id-dimensional lattice nearest-neighbor Laplacian: 

(A*)(x)= J2 (*(y)-*( x ))= E *(y) - 2dn*(x), (2.1) 

|y- x li=i |y- x li=i 
acting on £ 2 (Z nd ) : * G ^ 2 (Z nd ), x G Z nd . 

2.2.2. Assumptions on U. 

(I) Short-range interaction. The potential energy of inter-particle interaction 

U: (Z d ) n -> M 

is of the form 

U(x)= ^=(xi,...,x n )eZ nd (2.2) 

l<i<j<n 

where 3> : Z d x Z d -> K is a symmetric bounded nonnegative function which satisfies 

there exists ro G [0, oo) such that 

|xi -x 2 \ > r => $(xi,x 2 ) = 0. (2.3) 
We will call ro the "range" of interaction U. 

Some auxiliary results, concerning eigenvalue concentration (Wegner-type estimates), 
do not require the interaction U to have short range. Hence, it is convenient to formu- 
late the following, weaker assumption on U. 

(10) Bounded interaction. The potential energy of inter-particle interaction 

U: {Z d ) n -> M 

is of the form (|2.2[) . invariant by permutations of the particles x±, . . . , x n , and bounded. 

2.2.3. Assumptions on V . To describe our assumptions on the random potential 

V : Z d x Q -> R, 

we need some notations. Let Q C Z d be a bounded parallelepiped. We denote by 
£q(cj) the sample mean of the random field V over Q: 

i Q (u) = {V) Q = ±J2 V M, 

and by r) x (uj) the fluctuations relative to the sample mean 

r) x (u) = V(x,uj) -£q(w), 
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Let Fy be the common probability distribution function of the random potential V: 

F V (E) = F{V(0,lo) < E}. 

Now denote by 3V,q the sigma algebra generated by {r] x ,x G Q} and by ■ |3V,q) 
the conditional distribution function of £q given 3v,Q : 

Our assumption on V is as follows: 

(PI) External (random) potential energy. The random field V : 7L d x O — > R is 
i.i.d., nonnegative, and for any R > there exists a function z/r: R + — > R + with the 
following two properties: 

(i) for some b > 0, < i < 1 => j/ h (*) < C'R c "t b , 

(ii) (CM(i/)) For all Q C Z d with diam(Q) < i?, the conditional distribution function 

■ \$v,q) satisfies, for all i,s£l: 

eaasap\F € (t\Sv tQ )-F c {a\Sv t Q)\ < v R {\t - s\). (2.4) 

We also assumed that for all e > we have ¥{V(0,ui) < e} > 0, i.e., inf suppFy = 0. 

In Section [5j this assumption (PI) will be replaced by the following one: 

(P2) External (random) potential energy. The random field V : H l x £1 — > R is 
i.i.d., and the probability distribution function Fy satisfies 

supOFV(a + e) - F v (a)) < — ^ (2.5) 
aeR |lne|^ 

for C e (0, oo) and A > 24Nd. 
Note that this last condition depends on N. 

Remark 1. In the course of the scale induction, the bound (|2.5[) will be used in a 
situation where e = e~ L ,(3 = 1/2 and L > is a large integer. With such a value of 
e, (|2.5p takes the form 

sup(Fy(a + c- Lf< ) -Fy(a)) < Cr A . (2.6) 

Such a power-law estimate is sufficient for the purposes of the MSA: see p.4[) - (|3.5p . 
Assumption (PI) will give rise to much stronger bounds of the form e~ L 1 +°( lni ). 

2.3. The structure of the induction on the number of particles. Following the 
general approach of PQE] , we prove localization for systems with a number of particles 
bounded by some integer N < oo. The value of N > 2 may be arbitrary, but once it is 
chosen, it is fixed for the rest of the scaling analysis, and several important parameters 
used in our scheme depend upon N. 

As a result, we carry out a finite induction on the number of particles n varying 
from 1 (where the well-known results of the 1-particle localization theory can be used) 
to N. In some notations we emphasize the TV- dependence of parameters, while in other 
instances this dependence is omitted for simplicity. 

The inductive step from n — 1 to n requires information on n'-particle systems for 
all 1 < ml < n — 1. This can be briefly explained by the fact that the induction step 
involves the analysis of the decompositions of an n-particle system into two subsystems 
with n' and n" = n — n' particles, for any 1 < n' < n — 1. So it is assumed that all 
necessary bounds are available for n'-particle systems with any 1 < n' < n — 1. 
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Since our main method is (a multi-particle adaptation of) the MSA, the induction 
step from n — 1 to n involves also a scale induction, i.e., induction over a sequence of 
scales Lk, k > 0. The starting point of the scale induction requires key bounds of the 
MSA to be established at the initial scale i - 

2.4. Multiscale analysis. According to the general structure of the MSA, we work 
with lattice cubes 

c£°(u) = {x G 7L nd : |x - u| < L} . (2.7) 
The superscript, indicating the number n of particles, will be sometimes omitted, e.g., 
for n = 1. Define its internal boundary 

d-C^\u) = {v G 7L nd : dist (v, % nd \ c£°(u)) = l} , (2.8) 

its external boundary 

d+C^iu) = {v G Z nd \ c£°(u) = dist (v,C^(u)) = l} . (2.9) 

The cardinality of a cube C^ } (u) is |C^ l) (u)l = (2£ + l) nd - For simplicity, we will 
often use the simpler estimate |C^(u)| < (3L) nd . Let an n-particle cube 

c£°(u) - Ci(«i) x • • • x C L (u n ) c Z nd . 
We define its projection by 

nC ( ^(u) = C L ( Ul ) U • • • U C L {u n ) C 7L d . (2.10) 
More generally, for any nonempty subset J C {1, . . . , n} we define its j7-projection by 

njcM(u)= |Jc L K)cz d . (2.ii) 

We define the restriction of the Hamiltonian H^™) to a cube C}™ (u) by 



H (n) =H (n)| 



with Dirichlet boundary conditions on d + C}™ (u). 
We denote its spectrum by <r(H^ n ^ ) and its resolvent by 

G o£ o (u) (E) := ( H W )(u) - E) ~\ E G R \ <r{v%lj • (2-12) 
Its matrix elements G c („)^(x, y; E) are usually called the Green functions of the 
operator H^. 

A permutation r G <5„ acts on a configuration x = (x\, . . . ,x n ) G Z nd of n parti- 
cles by tx = (a; T -i(i), . . . ,ar T -i( n )). We define the symmetrized distance between two 
configurations x, y G Z nd by 

<fe(x,y)= min |x-ry|. (2.13) 

Definition 1 (^-distant). Let A > 1 be fixed. Two subsets A, B C Z nd , 1 < n < N 
are called ^-distant if 

d s (A,B) > ime. 

The multiscale analysis is based on a length scale {Lk}k>o which is chosen as follows. 
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Definition 2 (length scale). The length-scale {£fc}fc>o is a sequence of integers defined 
by the initial length-scale Lq > 2, and by the recurrence relation 

Lk+i = [L%], 

where 1 < a < 2 is some fixed number. In this paper, a — 3/2. 

This length scale {Lk}k>o is assumed to be chosen at the beginning of the multiscale 
analysis, except that in the course of the analysis it is often required that Lq be large 
enough. 

Definition 3 (E'-resonant). Let n > 1, (3 — 1/2, and E G K be given. A cube 
c£° (v) C 1 nd of size L > 2 is called ^-resonant (E-R) if 



dist 



E,a(H 



(n) ) 
C«(v)J 



<c~ Lf> . (2.14) 



Otherwise it is called E-non- resonant (i£-NR). 



The next definition depends on the parameter a > 1 which governs the length scale 
of our multiscale analysis. 

Definition 4 (^-completely non- resonant). Let EeKbe given. A cube c£°(v) C 
jnd Q £ g j ze x > 2 is called ^-completely non-resonant (J5-CNR) if it does not contain 
any E-R cube of size > L x / a . In particular, Cl(v) is itself £J-NR. 

Definition 5 {{E, m)-singular). Let m > and £elbe given. Let N > 1 be fixed. 
A cube C^ l) (u) C Z nd , 1 < n < N is called (E, m)-nonsingular ((£,m)-NS) if 

max G c („ )f Ju,v;E) < e~^ m ' L ^ L , (2.15) 
vea-cW(u) L 

where 7(771, L, n) = 7(m, L, n, N) is defined by 

7 (m, L, n) = m(l + £ -V8)JV-n+i > TO , ( 2 .16) 
Otherwise it is called (i?, m)- singular m)-S). 

2.5. Multi-particle localization at low energies. We prove our main result on 
multi-particle spectral localization under two different sets of assumptions: 

• (K), (I), and (PI), 

• (K), (I), and (P2), which is less restrictive than (PI). 

As we have already pointed out, assumption (PI), introduced in [3], leads to a simpler 
proof (compared to the scheme described in [3]) and to more optimal decay bounds for 
eigenfunctions, especially in finite cubes. The latter improvement does not manifest 
itself in the qualitative statement on pure point spectrum and asymptotic exponential 
decay of eigenfunctions. We discuss this issue in more detail in Section [5] 

Theorem 1 (multi-particle localization at low energies). Let HW(a;) be a random 
N -particle Hamiltonian of the form 

N 

HW(w) = - A + ^Vfo.wJ + U, 

where A, U satisfy (K), (I), respectively and let 

E° = inf o-QH(w)) P-a.s. (2.17) 
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(A) If the random potential V satisfies (PI) , or if it satisfies (P2) and is bounded from 
below, then there exists E* > E° such that with F '-probability one, 

(i) the spectrum ofH(uj) in (E° , E*] is nonempty and pure point, 

(ii) any eigenfunction \&i(x, u;) with eigenvalue Ei(uj) £ [E°,E*] is exponentially de- 
caying at infinity: there exist a non-random constant m > and a random con- 
stant Ci(uS) > such that 

|*i(x,«)| <C i (w)e- m W. (2.18) 

(B) If the random potential V satisfies (P2) and is unbounded from below: 

iaf{E e R : F V (E) > 0} = -oo, 

then, for any positive m, there exists E* (ra) S R such that (i) and (ii) are valid for 
energies E < E* (m) with the given m as rate of decay. 

Proof. See Sections I3H1 for (A) under (PI), and Section [5] for (A) under (P2) and for 
(B). □ 

2.6. Mult i- particle localization under weak interaction. To formulate our result 
on stability of multi-particle localization under perturbations by weak interactions, it 
is convenient to write the iV-particle Hamiltonian as follows: 

HM(u) =H<5H = - A + gV(uj) + hU, g,heR, (2.19) 

where V(w) = V(x,cj) = V{x\,uj) + • • • + V(xn, oj) and h is a parameter measuring 
the strength of the non-random interaction U satisfying (I) (boundedness and short 
range). For h — 0, one has a system of N non-interacting quantum particles subject 
to a common random potential V : Z d x £1 — s- R. 

Technically, we prove in this paper stability of the MPMSA bounds obtained in 
[B], for the principal - and most interesting - application is to one-dimensional multi- 
particle systems with arbitrarily small disorder amplitude \g\ > 0. As it is well-known, 
Anderson localization occurs in one-dimensional lattice systems with any non-constant 
i.i.d. random potential: no regularity of its marginal probability distribution is required. 
In particular, the i.i.d. potential may even take a finite number K > 2 of values. So, we 
prefer here a more general assumption (P2) (log-Holder continuity of the marginal dis- 
tribution). In this connection, we would like to emphasize that more singular potentials 
(e.g., those taking a finite number of values) cannot be treated by existing methods 
(neither MPMSA nor MPFMM), since these methods treat iV-particle Hamiltonians 
with N > 1 as iVcZ-dimensional lattice operators, and specifically one-dimensional an- 
alytic techniques (products of random matrices and Furstenberg theory) fail to apply 
to such models. 

Theorem 2 (multi-particle localization under weak interaction). Let g,h 6l and let 

hJ'h^aijVH + au 

be an N -particle hamiltonian in i 2 (1 Nd ), with A, U, and V satisfying (K), (I), and 
(P2), respectively. 

We assume that H^ } (o;) ; i.e., the Hamiltonian without inter-particle interaction, 
satisfies the following condition: there exists a nonempty energy interval I C M such 
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that, for all E G I , 

P{c£?( u ) is (E,m)-S) < L~ 2p , (2.20) 



forp>6Nd, m>CL 1/2 , 

and with Lq large enough. 

Then there exists h* > such that for any h G (—h*, h*) the Hamiltonian (w), 
with interaction of amplitude \h\, exhibits Anderson localization in the interval I, i.e., 
with ^-probability one, 

(i) the spectrum o/H^ (w) in I is pure point, 

(ii) the eigenfunctions \&j(x, w) relative to eigenvalues Ei(w) G / are exponentially 
decaying at infinity: 

|*i(x,w)|<Ci(w)e- m|x l m>0. 

Proof. See Section RTT1 □ 

The following statement relative to A-particle systems in Z will be derived from 
Theorem [5] since, for d — 1, we are able to establish the estimate (|2.20[) . 

Corollary 1 (TV-particle localization under weak interaction, d = 1). Let 

H 2M = -A + .gV(x,w) + /iU, g,heR 

be a random operator in £ 2 {1 N ) describing an N -particle quantum system in the one- 
dimensional lattice Z subject to the external random potential g\(uS) and to the inter- 
particle interaction hXJ , with A, U, and V : Zx(]->R satisfying (K), (I) ; and (P2), 
respectively. 

Then for any g =/= there exist h* > and a nonempty interval I C 1R such that 
Hg N ^(uj) exhibits complete Anderson localization in I for all h G (—h*,h*), i.e., with 
V '-probability one, 

(i) the spectrum o/H^(cj) in I is pure point, 

(ii) there is an orthogonal basis {i/ji(x, oj)}i £ N inl 2 (7j N ) consisting of eigenfunctions 
ipi(x,u)) relative to eigenvalues Ei(uj) G / with exponential decay at infinity: 

\Mx,u)\ <a(w)e- m|11 , m>0. 
Proof. See Section W% □ 

Taking into account that our approach is perturbative, it is clear that the value of 
the "mass" m — m(g, h, Fy) > 0, depending upon the PDF Fy of the random potential 
V(lu), for \h\ small enough is close to the upper Lyapunov exponent of the corresponding 
single-particle one-dimensional Anderson model. For bounded potentials and small \g\, 
the latter is known to be of order of 0(g 2 ) for energies E away from the center and 
edges of the interval [0, 4] giving the spectrum of the nonnegative lattice Laplacian — A 
in Z (cf. [7, 13, 14J and references therein). 

3. The TV-particle MSA scheme under (K), (I) and (PI) 

3.1. The general structure of the MPMSA. Recall again that we always assume 
N > 2 to be chosen and fixed. We consider only systems with no more than N particles. 
The proof is by induction on two parameters, the number of particles n = 1, . . . , N 
and the scale index k > 0. See section HOI 



MULTI-PARTICLE ANDERSON LOCALIZATION AT LOW ENERGIES 



9 



Assume that Lq > 8tq. Given 1 < n < N, the following property for some E* > E° 
will play an important role in our main strategy. 

(T>S.k,n, N). For any pair of L^-distant cubes C^^u) and C^(v) 

p{3£eJ: C£j(u), c£J(v) are (E,m)s\ < i^ 4 "'" (3.1) 
where m > 0, p > 6Nd, and / = [E°, E*] with E* > E°, are fixed. 

By Definition [U if for N fixed and 1 < n < N two cubes C^(u) and C^(v) 
are L^-distant then |u — v| > 11 NLk- The role of this lower bound on the distance 
between their centers will become clear later. We denote (DS.k,N) := (DS.fc, N, N). 

(DS.fc, A). For any pair of L^-distant cubes C^p(u) and (v) 

P{3E e I : cg?(u), Cj»(v) are (E,m)s} < (3.2) 

where m > 0, p > 6 Ad, and I = [E°, E*] with E* > E°, are fixed. 

Theorem 3. Let H^fw) = - A + Y,f =1 V{ Xj ,uj) + U, where A, U, V satisfy (K), 
(I), (PI), respectively. 

Then for any p > there exist E* = E* (p, A) > E° and m > such that (DS.fc, A) 
holds true for all k > provided Lq is large enough. 

Proof. The proof occupies Sections [3] and |3J □ 

Theorem 4. Let HW(u) = -A + £^Li Vfo.a;) + U, where A, U, F satisfy (K), 
(I), (PI), respectively. Suppose (DS.fc, A) ZioWs /or some _E* > and all k > 0. 
TTien wzi/i ^-probability one, 

(i) t/ie spectrum ofH^ N '{ui) in [E°,E*] is pure point, 

(ii) i/iere exists a non-random number m > suc/i t/iat eigenf unctions (x, w) o/ 
HV^a;) wrtfft eigenvalues Ei(u>) £ [E°,E*] are exponentially decaying at infinity 
with rate m: 

|*i(x,w)| <CiHe- m l x l. (3.3) 

Proof. The idea of this theorem goes back to the pioneering works [8, 10J and it can be 
proven essentially in the same way as in these papers, with minor adaptations. See, 
e.g., 0. □ 

It is readily seen that Theorem 2] combined with Theorem [3] implies Theorem [T] (A) 
if V satisfies (PI). It remains, therefore, to prove Theorem [3] 

3.2. Eigenvalue concentration bounds. Given n > 1 and Lq > 2, we introduce the 
following two properties (WSl.n) and (WS2.n) of random Hamiltonians H c („)^. 

(WSl.n). For any L > Lq, any cube C}™ (x), and any E eR, 

P{c£°(x) is not £-CNR} < L~ q . (3.4) 
(WS2.n). For any L > Lq and any pair of L-distant cubes (x) and C^(y), 

P^3E e R, neither c£°(x) nor C^ } (y) is £-CNr} < L~ q (3.5) 
for some q = q(n) large enough. 
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More precisely, for a fixed number of particles N, and any 1 < n < N , we always 
assume that the parameter p > in (DS.fc, n, N) or (DS.fc, N) and the exponent q in 
(WSl.n) or (WS2.n) satisfy 

p > 6Nd, q > 4 N p > 6 x 4 N Nd, (3.6) 

i.e., in terms of the parameter a (we set a — 3/2): 

p > 4Nda, q > i N p > 4 N+1 Nda. 

Lemma 1 ([!]). Under assumptions (K), (10), and (PI) or (P2), both properties 
(WSl.n) and (WS2.n) hold true for any positive integer n provided Lq is large enough. 

Lemma 2. Under assumptions (K), (10), and (P2), for any cube (u), L > 1, 
1 < n < N, and any s G (0, 1] 

p{||G cg0(u) (£)| >.^} <C|C L (u)| 2 x \ln S \~ A (3.7) 

where A > 247V d is as H2.5|> . which is part of assumption (P2). 

3.3. Initial scale estimates. Let N > 1 be fixed and let 1 < n < N. We observe 
that (DS.0,n, N) follows from the following stronger property: 

(S.O, n, N). For any x e Z nd , 

P ja£ e [E°, E*] such that C^(x) is (E,m)s\ < L~ 2piN ~ n (3.8) 

where m > 0, p > 6Nd, and / = [E°, E*] with E* > E°, are fixed. 

We denote (S.O, TV) := (S.O, N, N). 

(S.O, For any x e Z Nd , 

¥^3E e [E°,E*] such that C^(x) is (E,m)s\ < L a 2p (3.9) 

where m > 0, p > 6Nd, and / = [E°, E*] with E* > E°, are fixed. 

Indeed, (DS.O, n, A^) assesses the probability to have two cubes to be singular at the 
same energy. Such an event implies that the first cube (as well as the second one) is 
singular, and for one cube one can apply (S.0,n, N). So, for our purposes, it suffices 
to establish (S.O, n, N) for any 1 < n < N. 

The initial scale bound (|3.8p — hence also (DS.0,n, N) — will be proven with 
the help of a similar property of single-particle operators. Specifically, we need the 
following two results from the single-particle localization theory. 

Lemma 3 (Combes-Thomas estimate). Consider a Schrodinger operator 

H = -A + W(x) 

acting in £ 2 (A), A C 1? , v>\, with an arbitrar$\ potential W : A — > K. Suppose that 
• E eM is such that dist(S, cr(-ff)) > r\ with < r) < 1. 
Then for all x,y G A 

(H - E)- 1 (x,y) <-e-^\ x -y\i<- e -Th\ x -y\. (3.10) 
r, r) 

Proof. See e.g. [Ill Theorem 11.2] □ 



^This includes the cases of single- and multi-particle operators, that differ only by their potentials. 
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Lemma 4. Let HW(u) = -A + V(x, u) be a random single-particle lattice Schrd- 
dinger operator in £ 2 {lj d ). The random variables V(x,ui) are i.i.d., non constant and 
nonnegative. 

Then, for any C > and arbitrary large Lq > there exists c > such that the 
lowest eigenvalue Eq^^uj) of the restriction (u)( w ) °f H^(oj) to a cube Cl (u) 

with Dirichle^ boundary conditions satisfies 

P{4 1} M < 2CL- 1/2 } < e - c l^o(«)l 1/4 . (3.11) 

Proof. See [111 Theorem 11.4] and its proof. Lemma U follows actually from the study 
of Lifshitz tails, see [TTJ Lemma 6.4] and |1 51 Theorem 2.1.3]. □ 

In our earlier work [9], we have already inferred the initial scale MPMSA estimate 
for the two-particle lattice Anderson model from the results given by Lemma [3] and 
Lemma The adaptation to an arbitrary number of particles n is straightforward. 

Lemma 5. Let H( n '(cj) = —A + V(xi,u>) + • • • + V(x n ,u>) + U(x) be an n-particle 
random Schrddinger operator in ( 1 {fL nd ). Suppose 

(i) the interaction U is bounded and nonnegative^, 

(ii) the random external potential V{x,uj) is i.i.d. and nonnegative. 

Then, for any C > and arbitrary large Lq > there exists c > such that the lowest 
eigenvalue Eq(w) of the restriction H^ n ;L (w) ofH.^ n \uj) to a cube C^(u) C 

C Lq ( U ) 

satisfies 

P{4 n) H < 2CL- 1/2 } < e - c l^o(«)l 1/4 . (3.12) 

Proof. Since the interaction potential U is nonnegative, then by the min-max principle 
the lowest eigenvalue Eq(oj) of H^"^ is bounded from below by the lowest eigenvalue 

Eq U \uj) of the operator 

HW(u) = -A + V(xi,u) + --- + V(x n ,u). 
The latter operator can be rewritten as follows: 

n 

H (,l) (tj) = V l £ 2 (Z d) ® ■ • • ® l £ 2 (Z d) ®Hf ] ® V (Z d) ® • ■ • ® l f 2 (Z d) 

j = l * * 

J—1 times ri— J times 

where Hj — — -f V(xj,o;), j = 1, . . . , n. Hence the lowest eigenvalue Eq(cj) of 
HW(w) has the following form: 

n 
3=1 

nci ^ V>c» lrviTT-oc-f- ni rrnmrolno FT"* * All oirromrQlnno K7*^' 



where £Jq . is the lowest eigenvalue of H y - ' . All eigenvalues Eq • are nonnegative due 

to the nonneg, 
for any s > 0, 



to the nonnegativity of the operators with nonnegative external potential, so that 



F{Ei n \u)<s}<F{E$(uj)<s}. 



2 In the proof, one usually considers first an operator with Neumann boundary conditions, and then 
applies Dirichlet-Neumann bracketing. 

"^The assumption that U has short range, or is decaying at infinity, is not required. 
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Therefore, Lemma [5] follows from Lemma 0] applied to the single-particle Schrodinger 
operator H± . □ 

Remark 2. For the assertion of Lemma [5] to be useful, i.e. in order to ensure that the 
spectrum of random operators tv- n '(uj) in £ 2 (Z rad ) is nonempty in any interval [0,e], 
e > 0, one would need to require that the marginal probability density function Fy 
satisfies 

£* := M{E e M : F V {E) > 0} = 0. 
Formally, however, the assertion of the lemma remains valid even for £* > 0. In that 
case, owing to nonnegativity of the operators — one has > £* > 0, so that, 
trivially, 

p|^oH < = 0. 

Theorem 5 (initial scale estimate). Let N > 1 be fixed and HW(u;) = —A + 
V{x\ 1 lj) + • • • + V(xn, w) + U, where A, U, V satisfy (K), (I), (PI), respectively. 

Then for any p > 6Nd, for any Lq large enough and < C < oo, there exist 
E* > E° and m > CL > such that properties (S.0, n,N) and (DS.O, n, N) hold 
true for any 1 < n < N . 

Proof. Apply Lemma [3] and Lemma [S] □ 

1/2- • 1/2 

Remark 3. Clearly, the inequality m > CL implies m > CL k for all k > 0. 
Further, by definition of the function 7(m, L) — see (|2.16l) . for all L > and 1 < n < N, 

— 1/2 

"f(m, L, n) > m. Therefore, if one proves (DS.O, n, N) with m > CL 1 , then the lower 
bound on the decay exponent of Green functions, required at each subsequent step k 
of the scale induction, will be guaranteed. 

For future use we also give the next statement which treat the case of external 
random potentials which are bounded from below. 

Lemma 6 (V bounded from below). Let H> n ) (u) = —A + V(xi , w) H VV{x n ,uj) + 

U(x) be an n-particle random Schrodinger operator in l 2 {T, nd ). Suppose 

(i) the interaction U is bounded and denote 

|U = inf U(x), 

(ii) the random external potential V (x , uj) is i.i.d., satisfies (P2), and is bounded from 
below: there is £ * E K such that the probability distribution function Fy satisfies 

m£{E e M : F V (E) > 0} = £ *. 
Then, for any C > and arbitrary large Lq > there exists c > such that the lowest 
eigenvalue Eq(lo) o/H ( „) (w) satisfies 

^Lq ( U ) 

P [Eq(u) <n£*+^ + 2CL Q 1/2 } < e - c l Ci o(«)l 1/4 . (3.13) 

Here, H ( n) (ui) denotes the restriction o/H^™^(w) to the cube C^^(u) C II nd . 
®l ( u ) 

Proof. Let us set 

V(x,lu) = V{x,u)-£* 
U(x) = U(x) - p. 
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Since V(x,ui) > £* almost surely and U(x) > fj, we have 

V(x,uj) > and U(x) > 0. 

Therefore, 

n 

H<£ j (a,) = - A + £ n^.w) + U(x) 
t0 j=i 

= - A + 7(0:1,0;) H h V(x n ,u) +U(x) + n£* + (i 

L o 

so that the lowest eigenvalue Eq(lo) of the random Hamiltonian H ) n) (lu) satisfies 

E (u)) = E (u}) + n£* + [i, 
where Eq(lo) > is the lowest eigenvalue of n ) n) (lo) + U(x). Next 

^L ( U ) 

E (lo) < n£* +{i + 2CL~ 1/2 E q (lo) < 2CL~ 1/2 . 

Thus, 

P {Eq(uj) < n£* + ft + 2CL 1/2 } < P (e (lo) < 2CL~ 1/2 } 

< e -c|CL («)| 1/4 _ 

Since U(x) and V(x,lo) are nonnegative, Lemma |H] follows by Lemma [S] □ 

Now we also treat the case where the random potential V is unbounded from be- 
low; probably, the most important application is to Gaussian i.i.d. random fields. A 
particularity of this class of random potentials is that the "mass" m > here can be 
made arbitrarily large, in the energy band {E < E* < 0} with sufficiently large \E*\, 
depending upon the desired value of the decay exponent m. Also, the proof of the 
initial scale bound in such a case is much simpler: one can use essentially the same 
argument as in the case of "strong disorder", without appealing to the Lifshitz tails 
asymptotics. 

Lemma 7 (V unbounded from below). Let HW (lo) = — A.+V(xi,u)-\ VV(x n ,io) + 

U(x) be an n-particle random Schrddinger operator in l 2 (Z nd ). Suppose 

(i) the interaction U is bounded, 

(ii) the random external potential V(x,lo) is i.i.d., satisfies assumption (P2), and is 
unbounded from below: 

M{E £ R : F V (E) > 0} = -00. 

Let rj > 0, q > and Lq > 2 be given. Then there exists E* = E*(rj,q, Lq) < such 
that the lowest eigenvalue E^(lo) o/H^"^ (lo) satisfies the estimate 

P{4"V) <E*+T]} < LP- (3.14) 
Here, H^, (lo) denotes the restriction ofH^(u}) to (^(u) C Z nd . 
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Proof. First of all, note that, by nonnegativity of the kinetic energy operator and 
boundedness of the interaction U, 

-A + U> -C I 

with < C := ||U|| < oo. Therefore, for any s < and any fixed Lo G N, 
F \E { n) < ns\ < p{ min V(x, lj) < ns + c\ 

1 J l x£C<;»(u) J 

<p{ min V{x,u)) < s + Cn~ x \ 

< |77C ( r n) (u)l x F v {s + Cnr 1 ) ► 

where V(x; u>) = V(x\,u) + • • • + V{x n ,oS). Thus for any fixed integer Lq > 2, any 
rj > and q > there exists E* such that (|3.f4l) holds true. □ 

3.4. The scale induction step. It remains to derive (DS.fc+1, n, N) from (DS.fc, n' , N) 

for any 1 < n' < n, i.e., to prove the following statement: 

Theorem 6 (scale induction step). Let HW(u;) = — A + YljLi ^(^ji^- 1 ) + U, where 
A, U, V satisfy (K), (I), (PI), respectively. 

Then there exists < L* < oo such that, if Lq > L* , then, if for k > 0, property 
(DS.fc, n', N) holds true for any 1 < n' < n < N , then (DS.fc + 1, n, N) holds true. 

Proof. This theorem splits into Theorem [71 Theorem [51 and Theorem [51 The proofs 
of these theorems occupy Section 2) Their end will mark also the end of the proof of 
Theorem [TJ (A) under (PI). □ 

4. MULTISCALE INDUCTION UNDER (K), (I) AND (PI) 

In section l4"Jl we establish some useful geometrical facts valid for any n > 1. Starting 
from section l4~2l we perform the induction step n— 1 ~^ n, using scale induction in L^, 
assuming all necessary properties of n'-particle systems with any 1 < n' < n — 1 and 
in cubes of any size Lj, j > 0. Further, the induction step L^ L^+i is carried out 
separately for three types of pairs of singular n-particle cubes (see sections I4.2M. 41) . At 
this step, we are allowed also to assume property (DS.fc, n,N) of TV-particle systems 
in cubes of size Lk, and aim to prove (DS.fc + 1, n, N) in cubes of size Lfe+i. 

As indicated by the title of this section, we assume here the condition (PI) which 
gives rise to a simpler MPMS A scheme than in [5] . The key point is a more convenient 
eigenvalue concentration bound (|3.5p (see Lemma [TJ) established earlier by Chulaevsky 
[4J. However, this bound has not been proven so far for a more general class of i.i.d. po- 
tentials with log-Holder continuous marginal distribution. For this reason, we provide 
in Section [S] an adaptation of the original MPMSA scheme used in [B] (for strongly 
disordered systems) to the localization at low energies, under a weaker assumption 
(P2). 

4.1. Fully and partially interactive cubes. Consider the diagonal in Z nd , n > 1: 
ID={x=(x,...,x)£Z nd :x£Z d }. (4.1) 

Definition 6 (fully/partially interactive). An n-particle cube C^(u) C Z nd is 

(i) fully interactive (FI) if dist(u,D) < 2n(L + r Q ), 

(ii) partially interactive (PI) otherwise, i.e., if dist(u, D) > 2n(L + tq). 
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The following simple statement clarifies the notion of PI cube. 

Lemma 8. If dist(u, B>) > 2n(L + ro), then there exists a subset J C {1, . . . , n} with 
1 < card J < n — 1 such that 

max \uj 1 — uj 2 \ > tq. (4-2) 

hi j 

Consequently, if a cube C^(u) is PI, then there exists a subset J C {I, . . . , n} with 
1 < card J < n — 1 such that 

njci n) (u)nnj,c { l l \u) = 0. 

In other words, setting ri = \J\ > 1, n" = \J C \ > 1, u' = Uj £ Z dn ' , and 

^present (u) as a cartesian 

cW(u) = ci"''(u')xCf)(u") (4.3) 



u" =ujc e 7h dn , one can represent C^(u) as a cartesian product 



with disjoint projections 

njC { l l \u) = nc^\u')cz d , 
n JC &l l \u) = nc ( [ l "\u") c z d . 

For a given PI cube, there can exist several decompositions of the form (|4.3[) . We 
will assume that such a decomposition is associated with each PI cube and call it the 
"canonical" decomposition. 

Now we turn to geometrical properties of FI cubes. 

Lemma 9. Let x = (x x , . . . , x n ) £ Z d x • • • x Z d ~ Z nd . Then 

dist(x, D) < max \x, - Xj\ = diam 77x < 2 dist(x, ED). (4.4) 

i,j 

Proof. By dehnition of D, 

dist(x, H5) = min max \xj — u\ 

< min max \ xj — Xi \ (each Xi IS cL particular u G Z d ) 

i j 

< max \xj — Xi\ 

id 

= diam{a;i, . . . , x n } — diam Tlx, 

proving the first inequality. On the other hand, if u = (it, . . . , u) G ES C Z nd is the 
closest point to x in D, then, for any 1 < i,j < n, 

\xi -xj\ < \xi -u\ + \u -Xj\ < 2 dist(x, u) = 2 dist(x, ED). □ 

Lemma 10. Let n > 1, L > 8ro and consider two n-particle fully interactive cubes 
C ( [ l) (u) and C^ } (v) with \x - y| > 11 nL. Then 

nc ( £ ] (u) n nc { £ ] (v) = 0. (4.5) 



Proof. If for some R > 0, 



R< |x-y| = max |^ - % | 

l<J<n 
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then there exists 1 < jo < n such that \xj — Uj \ > R- Since both cubes are fully 
interactive, dist(x, D) < 2n(L + ro), dist(y, D) < 2n(L + ro), so we can use (I4.4j) for 
the centers x, y and write: 

\xj — Xi\ < 2 x 2n(L + ro), 
\yj -yj\<2x 2n(L + r a ). 
By triangle inequality, for any 1 < i, j < n and i? > 11 > 10 nL + 8 nro, we have 

\ x i ~ Uj \ ^ \ x j ~ Vjo\ ~ \ x jo ~ x i \ ~ \Vjo ~ Vj\ 

> 10nL + 8nr Q - 8n(L + r ) = 2nL. 
Therefore, for any 1 < i,j ' < n, 

mmdist(C L (xi), C L (yj)) > min \xi - yA - 2L > 2(n - 1)L > 0. 

This means that 

iici n) (x) n nc£\y) = |J (C L ( Xl ) n C L (y 3 )) = 0. □ 

i,3 

It is clear from the proof that the condition |x — y| > 11 nL could have been replaced 
by a more explicit, and sharper, inequality |x — y > 10 + 8nro, but we chose the 
first inequality for simplicity. 

Given 1 < n < N, and assuming (DS.fe, n', N) for any 1 < n 1 < n, we will prove 
(DS.fc + 1, n, N), separately for the following three types of pairs of cubes: 

(I) x (x) and (y) are both Pi-cubes. 

(II) C^^x) and C^ +1 (y) are both Fl-cubes. 

(Ill) One of the cubes ^x) or x (y) is PI, while the other is FI. 
In the rest of the paper we will denote by / the interval [E°, E*]. 

4.2. Pairs of partially interactive cubes. We consider here a pair of Lfc + i-distant 
partially interactive cubes C^ +i (x) and C^ +i (y) and prove that (DS.fc+1, n, N) hold 
for such a pair. We are allowed to assume property (DS.fc, n', N) for any 1 < n' < n 
and every pair of L^-distant cubes C^™ •'(u) and C^™ '(v). 

Let C^'^^u) be a Pi-cube with the canonical decomposition (|4.3[) . We also write 

x = (x',x") for any point x £ C}™ ^u), in the same way as (u',u"). In fact owing 
to the symmetry of the operator U, we can simply take without loss of generality 
J = {l,...,n'} and J c = {n' + l,...,n} so that the corresponding Hamiltonian 
H_,(„) , , is written in the form: 

C> (u) 



Hi, , .*(x) = (-A*)(x) + [U(x') + V(x» +U(x") + V(x",w)]*(x) (4.6) 



or, in compact form 

Here H , and H^' 1 . L are the Hamiltonians acting on x' and x", respec- 

C^' 1 > (u') 1 (u") 6 

tively, while I is the identity operator on the complementary variable. We denote by 
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q(" )(u',v';.E) and G^ n >(u",v";E) the corresponding Green functions, respectively. 
Let 

{( A ^i)} i=1 ,...,|C^ i( u')| and UM>^h=i,-,\C^ x (u»)\ 

be the eigenvalues and corresponding eigenfunctions of H^"^,, and H ( -™J, ) , 

c ^+i (u,) Ci 1+i (u " ) 
respectively. Then we can choose eigenvectors VPy of H ( n) as tensor products 

* y (x)=^(x')<8)^(x") (4.7) 
while corresponding eigenvalues 

/••'., A, • //,. (4.8) 

The eigenvectors of finite volume Hamiltonians appearing in subsequent arguments and 
calculations will be assumed normalized. Introduce the following 

Definition 7 (m-tunnelling). Let 1 < n < N, k > 1, I = [E°,E*] and m > be fixed. 

(i) A cube (u) is called m-tunnelling (m-T) if there exist E £ I and two L^-i- 
distant cubes (vi), C^" ) _ i (v 2 ) C C^(u) which are both (E,m)-S. 

(ii) A cube C^(u) is called (m,n' ,n")- partially tunnelling ((m,n',n")-PT) if it has 

a factorization C^'V) x C^'^u") in which either C^'V) or c z™'' V") is 
m-tunnelling. Otherwise it is called (m, n' , n")-NPT. 

(iii) A cube (u) is called m-partially tunnelling (m-PT) if it is (m, ?(', n")-PT for 
some n', n" > 1 with n' + n" = N. Otherwise, it is called m-NPT. 

The following statement is borrowed from [5] . 
Lemma 11. Let C^(u) be a PI cube of the form 

c£>(u) = c£V)xc£V) 

where u = (u', u") with u' = . . . , u n >) G Z™ ' d and u" = (v+i, ■ ■ • , u n ) £ Z n d . 
We assume that 

(i) \uj 1 — Uj 2 1 > 2Lfe + rg for all I < ji < n' , n' + 1 < j% < n. 

(ii) C^(u) is m-NPT for some m > 0. 

(iii) C^(u) ™ E-CNR for some Eel. 
Then, 



max 
i<z<|ci';'>| 

v"ea^c ( "" ) (u" 



max 

i<j<|c4 n fc ''V')l 
v'ea-c("''(u') 



G ( ^""- I (u", v"; E — \i) <e~ m ' Lk (4.9) 
G^'Vy^-Mi) <e~ m ' Lfc . (4.10) 



where 



m> = 7 (m, L fc , n - 1) (l - - InL^^j (4.11) 

and where and {/J-j} are the eigenvalues of H A) and H^" (n ,,j , respectively. 

c i™, ("') c i" ( u ") 
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Proof. See Lemma 3]. □ 

Lemma 12. A PI cube C^(u) which satisfies assumptions (i)-(ii)-(iii) of Lemma \ 1 1\ 
is(E,m)-NS. 

Proof. Using Definition (12 . 16(1 which gives the explicit form of "f(rn, Lk,N) t it is easy 
to see that m! > j(m, Lk,n), hence inequalities (|4.9[) and (|4.10l) of Lemma [TT] imply 



l<i<|C<™'>] 

Ved-c^piu") 



i<j<|c£™ »(u")| 
v'e8-ci"' ) (u') 



G( n ">(u",v";f?-Ai) < e - <{m ' Lk ' n)Lk , (4.12) 



Then C^(u) = Cj^(u') x C^" ^(u") is (E, m)-NS, using the same representations 
G W (u, v; E) = J2 ^(u')^(v')G( w ")(u", v"; E - A,), 

i 

GW (u, v; £) = £ Vi(u>i (v")G("')(u', V; E - H ) 

3 

as in [6, Appendix B]. □ 

Lemma 13. Let 2 < n < N. We assume property (DS.fc, n 1 ,N) for any 1 < n' < n. 
Let C}™ (y) be a PI cube and let 

C £ 1 (y) = c £! 1 (y')xc(:2(y") 

be a factorization, y = (y', y"), y' = (y 1} . . . , y„/) G y" = (y n /+i, . . . , 2/at) G Z"" d 

toit/i n — n' + n", n', n" > 1 and 

min |y; -%| > 2L fc + r . (4.14) 

l<i<n' 
n'+l<j<iV 



P{C^ +i (y) z S m-PT) < \liI^~\ (4.15) 
Proof. By definition we have 



Ty := 



{Cg +i (y) is m-PT} C {either C^y') or C^fy") is m-l) . 



Since n' < n we can use property (DS.fc, n', N). Thus by combining this property with 
an upper bound of the number of possible pairs of centres vi, V2 G Gj™ ^(y'), we get: 

P{ci";| 1 (y / ) is m-T} < ^LUt x L~^ N - n '. (4.16) 
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Similarly for P{C^™ \ (y") is m-T}, using (DS.fc, n" , N) which is known since n" < n. 
Therefore, 

P{T y } < P{either C^y') or C^y") is m-PT} 

< C(n,N,d)L k+1 

, ^/ Ar j\ T -2p A 4 JV -"+2(n-l)d 

< C{n 1 N 1 d)L k+1 = 

< C{n, N, d)L k * x (a = 3/2) 

< \^~" (4-17) 

for large enough Lo, since p > ActNd = 6Nd and g > 4rp by our assumption. □ 

Theorem 7. There exists L\ > smc/i f/iaf if Lo > L* and if for k > and 1 < n < TV, 
(DS.fc, n' , N) holds true for any 1 <n' < n, then (DS.fc + 1, n, N) holds true 

(i) for any pair of L^+i- distant PI cubes C^ +i (x) and i (y). 

Proof. Let C^™; t (x) and i (y) be two L^+i-distant Pi-cubes. We consider the 
events: 

B k+1 = {3EeI: Cg +i (x) and c£J +i (y) are (25,m)-S}, 
R = {3£ e I : neither (x) nor (y) is £-CNR}, 

T x = {C^ +i (x) is m-PT}, 
Ty = {c£ ) +i (y) is m-PT}. 

If u G B k+1 \ R, then V£ G J, c£J (x) or cg +i (y) is £-CNR. If c£J +1 (y) is £7- 
CNR, then it must be m-PT: otherwise it would have been (E,m)-NS by Lemma [T2l 
Similarly, if C^J (x) is S-CNR, then it must be m-PT. This implies that 

B fe+1 C R U T x U T y . 

Therefore, 

P{B fc+1 } <P{R}+P{T x } + P{T y } 
< L~l, + 2P{T X } 



-2p4 J 



J k+1 

where we used (WS2.rc) to estimate P{R} and (g7T7|) to estimate P{T X }. □ 

Given an n-particle cube C^ +i (u) and E G M, we denote by Mnd(C^"' I +i (u), _E) 
the maximal number of pairwise L^-distant, partially interactive (E, m)-singular cubes 
C^(uW) C C^ +i (u). Further, set 

M ND (c£> (u), 7) := sup M ND (C^ (u), 25). (4.18) 
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Similarly, we denote by Md(C}™ l (u), E) the maximal number of pairwise L^-distant, 
fully interactive {E, m)-singular cubes C^(u^) C C^™; ^u) and set 

M D (cg (u),/) := sup M D (Cg +1 (u),£). (4.19) 

fc+i B£/ &+i 

Lemma 14. W^it/i the above notations, 

p{M ND (cg +i (u),J) > 2} < (V + ^ 4P4 " _ ") • (4-20) 

Proof. The number of possible pairs of centres {u^, u^ 2 '} of L^-distant cubes 

c{!(uW), Cg(u^))ccg +i (u) 
which are PI and (_E, m)-singular is bounded by ^-Lijl]^, then 

P{A/ ND (C^ +I (u),7) > 2} < — LlH x P{B fe } 

with ¥{B k } < L~ q + 2P{T X } and P{T X } < ±L~ ipiN ~ n by (HTTl) . □ 

4.3. Pairs of fully interactive cubes. Our aim now is to prove (DS.fc + l,n,N) 
for a pair of Lk+i -distant fully interactive cubes C^ +i (x) and 1 (y). We always 
assume property (DS.fc, n, N) for any pair of ife-distant cubes. 

The main result of this section is Theorem [HJ We need some preliminary results. 

Lemma 15. Given k > 0, assume that property (DS.fc, n,N) holds true for all pairs 
of Lk- distant fully interactive cubes. Then for any I > 1 

p{m b (C { £ +i (u),I) > 2^} < C(n,N,d,£)Ll edna Lf ep4N ~\ (4.21) 

Proof. Suppose there exist 2£ pairwise ifc-distant, fully interactive cubes C^(u^^) C 
(u), 1 <j < 21. Then, by LemmaES for any pair C^u " - ^). C^u' 2 ' 1 ), the 

corresponding random Hamiltonians Hs', , and H*-™!, , are independent, 

i b c^V 2 - 1 ') c^(u(=«)W p 

and so are their spectra and their Green functions. Moreover, the pairs of operators 



H 



i = 1, . . . , £, form an independent family. Each operator H (w) (w), i = 1, . . . , 21 is 

indeed measurable with respect to the er-algebra Bi generated by the random variables 
{V(x,w) : x e iTC^^uW)}. By LemmaCTOl the projections UC^uW), i = 1, . . . , 2£ 
are pairwise disjoint, so that all sigma-algebras Bi, i — 1,...,2£ are independent. 
Thus any collection of events relative to the pairs (H ,L (a/), H L , (w)|, 

i = 1, . . . , £, also form an independent family. For i — 1, . . . , I we consider the events: 

A, = {3£e/: C^V 3 *- 1 *) and (u (M+2) ) are (£,m)-S} . 
Then by assumption (DS.fc, n, iV), we have, for i = 1, . . . ,£, 

P{A t }<L- 2p4N ~'\ (4.22) 



MULTI-PARTICLE ANDERSON LOCALIZATION AT LOW ENERGIES 



21 



and, by independence of events Ai, . . . , At,, 

p { n M = ri p (^) ^ ( L i 2p ) liN ~ n ■ ( 4 - 23 ) 



Ki<i 



To complete the proof, note that the total number of different families of 21 cubes 
C^" 3 (uW) C C^ +i (u), 1 < j < 21, is bounded by 

2£ 

C£ } +I (u) <C(n,N,£,d)Lf dna . □ 



(2£)\ 

Lemma 16 ([6, Lemma 12]). Fix an odd positive integer J and suppose further that 

(i) C^ +1 (u) is E-CNR, 

(ii) M^(C^ +i (u) : E)+M D (C^ k+i (u),E) < J. 
Then for Lq large enough, 



max 

vea-c'™' (u) 



G^„, ,.Ju,v;£) 



cy> (ay 

L k+1 v ' 



< e" m ' L 



where, for some C G (0, +oo), 

m' > 7 (m,L fc ,7i)(l-CL^ 1/2 ). (4.24) 

Using again the definition (|2.16[) of the quantity j(m, L,n), one can easily derive 
from Lemma [TBI the following bound. 

Lemma 17. Under assumptions and with notations of Lemma \lb\ (u) is (E,m)- 

NS. 

Proof. Corollary of Lemma [JJ)J Note that (|4.24p implies 

to' > 7(to, Lfc+i, n) (4-25) 
provided that Lq is large enough. □ 

Theorem 8. There exists > such that if Lq > Ljj; and if /or fc > and 1 < n < N , 
(DS.k,n, N) holds true, then (DS.fe + 1, n, N) holds true 

(ii) for any pair of L^+i- distant FI cubes i (x) and (y). 

Proof. Consider a pair of Lfc + i-distant cubes C^™; ^x), Cjj"; i (y) and set 

B fc+1 = [3E G I : c£> +1 (x) and C^ +i (y) are (E,m)s} , 

£ = [3E G I : neither C^ +i (x) nor C^ +i (y) is £-CNR,j , 

S X = [3EEI: M ND (C^ +i (x),i?) + M D (cgJ +i (x),B) > J+l}, 

S y = G J : M N D(C^ +i ( y ), 25) + M D (C^ +i (y), £) > J + l} . 

Let lu G Bfe + i, assume that lu £ £ U S x , then VE 1 G / either (x) or C^ n) (y) 

is E-CNR and M ND (C^ (x),£) + M D (C^ x (x),E) < J. Thus the cube (x) 
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cannot be _E-CNR: indeed, by Lemma [T71 it would be [E, m)-NS. So the cube 1 (y) 
which is 75-CNR and (E,m)-S. This implies again by Lemma [T71 that 

M NT3 (CtlM E) + M D (cg> +1 (y), E) > J + 1. 
Therefore w 6 S y , so that 

B fc+ i C SUS x US y . 

Hence. 



'{B fc+1 } <P{S} + P{S x } + P{S y } 

- -q 
J k+l 
— 4p 
J k+1 

Now let us estimate P{S X }. Set J = 21 + 1, with i — 2 then the inequality 



<L j ^ 1 + 2P{S x } 
<L-^ + 2P{S x }. 



M ND (C£ 3 +1 (x), £) + M D (C^ +i (x), £) > 2* + 2 

implies that either M ND (C ( ^ +i (x), E) > 2, or AfolC^ (x), E) > 21. Therefore, by 
Lemma [TJ] and Lemma [T51 

P{S X } <P{3££/: M m (C^ +i {x),E) > 2} 
+ P[3E 6 I : MdCC^^x),^) > 2^} 

< C(n,N,d) \ L k+ £ +L k+ \ I 

< C(n, /V, d) (£ fe 4 + L k+ \ j 

< C(n, iV, d) (L fc+ 1 + L fe - j [1 = 2, a = 3/2) 



< ~.E 



-2pi r 



4 fc+i 

where we used g > 4 w p and p > AaNd = 6Nd. □ 

4.4. Mixed pairs of cubes. Now it remains only to derive (DS.fc + l,n,N) in case 
(III), i.e., for pairs of n-particle cubes where one of the cube is partially interactive 
while the other is fully interactive. 

Theorem 9. There exists L%, > such that if Lq > L$ and if for k > and 1 < n < N, 
(DS.fc, n',N) holds true for any 1 < n' < n, then (DS.fc + l,n, N) holds true 

(iii) for any pair of Lk+i- distant cubes i (x), C^ +i (y) where one is partially 

interactive while the other is fully interactive. 



Proof. Consider a pair of Lfc + i-distant TV-particle cubes (x), (y) and sup- 

fc+1 



pose thatC^ n) (x)is partially interactive while C]™ (y) is fully interactive. Introduce 
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the events 

Bfc+i - {3E e I : C^ +i (x) and cg +i are (E,m)s} , 

S = {3E e I : neither C^ +i (x) nor C^ +i (y) is ^-CNr| , 
Tx = {c^ +i (x) is m-l), 

Sy = {BE e I : M m (C% +i (y\E) + M D (c£> +i (y),25) > J + l} . 

Let w G B fc+ i \(EU T x ), then for all E e 7 either x (x) or C^J (y) is £-CNR 

and (x) is m-NT. The cube (x) cannot be 7;-CNR. Indeed, by Lemma [TJ 

it would have been (75,m)-NS. Thus the cube ^y) is £-CNR, so by Lemma H71 

M ND (C^ ) +i (x),7) + A7 D (C^ ) +i (x),7) > J+l: otherwise C^ } +i (y) would be (75, m)- 
NS. Therefore ui e S y . Consequently, 

B fc+1 cSUT x US y . 

Recall that the probabilities P{T X } and P{S y } have already been estimated in Sections 
14.21 and 14.31 under the condition that (DS.fc, n', N) is known for any 1 < n' < n. Now, 
we obtain 

P{B fe+1 } <P{£} + P{T x } + P{S y } 

4.5. Proof of Theorem [T] (A) under (PI). Assuming property (DS.fc, n', N) for 
arr y 1 < n' < n we have proven (DS.fc + 1, n, N) 

(i) for any pair of Lfe + i-distant PI cubes (Theorem [7]) , 

(ii) for any pair of Lfc + i-distant FI cubes (Theorem [5]), 
(hi) for any mixed pair of 7/j+i-distant cubes (Theorem [9]) , 

i.e., we have proven (DS.fc + l,n, N) for any pair of Lfc+i-distant cubes: Theorem H2 
is proven. Since the initial scale estimate (for Lq) has been established in Theorem [5] 
both theorems imply Theorem [3] by induction in fc, hence assertion (A) of Theorem [TJ 

5. Proof of localization under (K), (I) and (P2) 

Now we will establish multi-particle Anderson localization at low energies under 
weaker assumptions on the external random potentials. Namely, we will not use the 
assumption (PI), the validity of which has been established so far for a more restricted 
class of random fields than (P2). 

To this end, we will provide an adaptation of the original scheme proposed in [BJ. 
However, this results in more complicated geometrical arguments and less optimal 
bounds on the decay of eigenfunctions. In particular, an eigenvalue concentration 
bound similar to Lemma Q] will not be established for all sufficiently distant pairs of 
cubes, but only for those satisfying an additional condition ("separable pairs"). This 
difficulty was encountered in [6j and in [T|. In the latter work, Aizenman and Warzel 
explicitly analyzed the nature of this technical problem and explained why they had 
to replace norm-distance estimates of the eigenfunction decay by those using the so- 
called Hausdorff distance; see the discussion in [2 §1.3]. A solution to this problem 
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has been proposed in [1], where the bound of the form (WS2.n) was proven under the 
assumption (PI), resulting in a streamlined description of the decay of Green functions 
and of eigenf unction. 

On the other hand, an adaptation of the multi-particle MSA scheme proposed in 
[5] to the proof of localization at low energies is fairly straightforward. Indeed, in the 
framework of the multi-particle MSA, like in the case of single-particle models, the role 
of the assumptions such as 

• strong disorder: \g\ 3> 1, or 

• low energies: E e [E°, E° + rj\, < T) <C 1. 

is to provide a mechanism which establishes the initial scale estimate; in our case, one 

— 1/2 

has to prove (DS.O, N) with a value of the "mass" to > CL and sufficiently large 
Lq. Apart from the fact that, unlike the case of strong disorder considered in [3], the 
value of the "mass" cannot be made arbitrarily large (to > m(g) — > 00 as \g\ — > 00), 
the rest of the MPMSA procedure can be performed exactly in the same way as for 

\g\ » i- 

For this reason, we only have to replace [3 Theorem 3] with Theorem [TUJ This 
theorem, combined with the arguments and results of [5] , will prove Theorem [T] under 
assumption (P2). 

For application to Hamiltonians with unbounded from below random external po- 
tentials, we give the following statement, namely the Combes-Thomas estimate for 
large spectral gap. 

Lemma 18 (Combes-Thomas for large gap). Consider a Schrddinger operator 

H = -A + W{x) 

acting in £ 2 (A), A C Z" ', v>\, with an arbitrary potential W : A — > M.. 
If dmt(E, (j(H) > 77 > 2, then for any n,m £ Z w , 

\{H - E)- l {n,m)\ < e~* . (5.1) 

Proof. The proof is obtained by a modification of that of Lemma [3] see the proof of 
[TTI Theorem 11.2]. □ 

Theorem 10 (initial length scale). Let HW(w) = - A + J2f=i v ( x j^) + U ; where 
A, U, V satisfy (K), (I), (P2), respectively. Let 

£* := mi{E € R : Fv{E) > 0} 

where F v is defined by F V (X) = P{V(0,u) < A}. 

(A) If £* > — 00, then for any p > 0, C > and arbitrary large Lq > there exists 
E* > E° such that properties (S.0, N) and (DS.O, N) hold true in I = [E°,E*} 
for some m > CL ^ 2 > 0. 

(B) If S* = —00, then for any p > 0, m > and Lq > 1 there exists E* > —00 such 
that properties (S.0, N) and (DS.O, AT) hold true in I = (— 00, E*] with the chosen 
value of m > 0. 

Proof. Taking into account the Combes-Thomas estimate (Lemma [3]) , assertion (A) 
follows directly from Lemma [B] in the same way as Lemma [5] in section f3 .31 

For (B), Lemma[7]and assumption inf{A e M : -FV(A) > 0} = — 00 imply that Ve > 
and V77 > 0, 3E*{ri) < such that 

P{£oM < E* +ry} < e. 
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On the other hand E (uj) > E* + r\ implies that for all E < E* 

dist(£,CT(if c (iv, (x) )) > T]. 

This allows to apply the above version of the Combes-Thomas estimate (Lemma [T5| 
with 77 > 2 and even arbitrarily large, that proves (S.0,7V), hence (B). □ 

6. Stability under weak perturbations 

Now we turn to the proof of multi-particle localization for systems with weak inter- 
action. 

Like in the case of disordered systems at low energies, the principal adaptation of 
the multi-particle multiscale analysis developed in [6] to weakly interacting systems, 
without the assumption of strong disorder (| g| ^> 1) lies in the verification of the initial 
scale estimate. 

Since the notions of Z?-resonance and (E, m)-singularity have been defined with 
respect to a given random ensemble of Hamiltonian, with fixed parameters, it is con- 
venient to change now the notations and to introduce explicitly the parameters g, h of 
H^(w). We will say that 

• a cube C^^u) is {E,g, /i)-resonant, if it is E'-resonant with respect to H^(w); 

• a cube C^^u) is (E, m, g, /i)-singular, if it is (E, m)-singular with respect to H^(w). 

6.1. Proof of Theorem [2j First of all we prove a lemma: 

Lemma 19. Suppose that the Hamiltonians Hg N J(uj) (without inter-particle interac- 
tion) fulfill the following condition: for sufficiently large Lq, 

p{3£eJ: C^(u) is (E,m,g,0)-S\ < L- 2p , (6.1) 

p > 6Nd, m > CL Q 1/2 , 

and that, consequently, the MPMSA allows to derive by induction the estimates 

p|3£6 7:C^(x) and Cf k \y) are (E, m, g, 0)~sj < L~ 2p , (6.2) 

for all k G N and any pair of L^- distant cubes Ci t (x), CL fc (y). 

Then there exists h* > such that for all h G (~h*,h*) the Hamiltonian Hg N ^(uj), 
with interaction of amplitude \h\, satisfies similar bounds 

p{3EeJ:Cg?(u) is(E,m',g,h)s} <L 2p ', (6.3) 

p' > 6Nd, ml > CL~ 1/2 , 

and consequently 

P e / : C^(x) and C^(y) are (E,m',g,h)-S\ < L~ 2p> . (6.4) 

Proof. Let us prove the assertion concerning the estimate (|6.3p . For simplicity, set 
m* = CL 1 ^ 2 with the same C, m, Lq as in (|6.ip . and 

G c - (u) , 9 ,J^) = (H^ )(u) ^-^r 1 , 5 ,^M. 
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V u; id y^, iu,y,vj-Lyu lll 

-1/8 



By definition, a cube C^ o ; (u) is (E, m, g, 0)-NS iff 

G c m (u)i} , (u,y;£) < e - m (i+L -v- )i0j (fi 5) 



where, by hypothesis, to > to*. Therefore, there exists e > such that 



max 
y69-c<"» ( u) 



1/Ss 



G c ,» )(u)9!0 (u,y;£) < e - m '^ +L o >*<> - e . (6.6) 



Since, by assumption, p > 6Nd, there exist 6Nd < p' < p and <5 > such that 
Lq 2p — 5 > Lq 2p . With such values p' and S, inequality (|6.ip with p > 6iVd implies 

P{3£ e I : Cj» (u) is (E, m, g, 0)-S} < L^' - <5. (6.7) 
Next, it follows from the second resolvent identity that 

ll G cr(u), ff .o(^) - G c<«W^)H < 1*1 l|U|| • l|G 0? o (B)iF|0 (£OII • HG cr(u) ^(^)l|. 

(6.8) 

By Lemma [21 applied to Hamiltoniant0 H^^ )( , ^ n and H^J]^ H u , for any <5 > 
there is -B(<5) € (0, +oo) such that 



n\\G ciViu)gfi (E)\\>B(5)}< S -, 
H\\G ciV{ ^ h (E)\\>B(S)}< 5 -. 

Therefore, 

P {HG <3(u)ff , (^)- G^^C^H > N||U||B 2 (5)} 

< P {H G ci«H 9 ,o^)ll >^)}+P{l|G c (-) (u)9 , h (^)ll >B(J)} 

„ * 5 
< 2- = - . 

" 4 2 



Set 



2||U||(B(d))2 >0 ' 



We see that if \h\ < h*, then 

\h\x\\U\\x(B{6)f<^. 

Hence, 

P {l|Gc<«»(u) s ,o-Gc Lo( u)^||>|}<2j. (6.9) 
Combining (|6.6p . (|6.8p . and (|6.9p . we obtain 
P{C^ } (u) is (£,m', 5 »-S} 

<P{cg?(u) is (B,m, ff ,O)-S}+P{||G c w (u)fl)O (£0-G c w (tt)flih (^)|| > |} 

<(L -^-5) + |<Lo V - □ 

Now Theorem [5] follows from Lemma [TH] combined with the assertion of Theorem 2] 
in which the interval [E°,E*] is replaced by /. 



^Note that the amplitude of the interaction has no impact on the bound 113.71 1. 
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6.2. Proof of Corollary [TJ In this subsection, d = 1. In order to apply Theorem [5] 
to weakly disordered A-particle one-dimensional lattice Anderson models, recall that 
one can prove sufficiently strong initial scale estimates for the Hamiltonian 

H{uj) = -A + gV{uj) 

with a non-constant i.i.d. random potential V: Z x 57 — > K (see Namely, for 

any g ^ and any interval /Cl, for some a = a(Fy,g) > and m = m(Fv,g) > 0, 
for all Lq large enough 

F{3E e I : C Lo (u) and C Lf » are (E,m)-S} < e~ aL °. (6.10) 

Actual proofs, which are non-perturbative in one dimension, allow to derive directly 
an exponential upper bound on eigenfunction correlators. Let 



then 



E 



T c („)(*, <e~ 3Q l^l. (6.11) 



It remains to prove that for the A-particle random Hamiltonian without interaction 
H {N \uj) in the one-dimensional lattice Z one has: 



P{cg?(u) ^ (E,m)-S}<L- 



■2p 



To this end we first give a preliminary argument dealing with some representation of 
the Green functions. Namely, 



G c ,» )(u) (u,y,i?)- ^_ _________ 

t, — 13 



,y ' ' ./^ E -(E jl +--- + E j 

iph ( U i)^h ivi) ■ ■ ■ i>jN ( u N)^j N (yw) 



E-(E h +--- + E jN ) 



Set 



ma * l G cf»(u)( u 'y^)l = \ G c ["\ u) ( u ,v;E)\ 

ye3-C<™'(u) L ° v ' L » v ' 



so that |u — v| = Lo, then there exists j £ {1, . . . ,N} such that \uj — Vj\ — Lq and 
since all eigenfunctions are normalised, we get: 



( u )\ T c L „(u 3 )(uj,Vj\u) 



\G n(N)( Ju,v;E)\ < 



dist(_,a(H™ )) 

sNl o t c Lo (u j )(u 3 ,v 3 ;uj) 
dist(_ >£ r(Hgjt, )) 
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Now, we have that: 



' 1 3E G I : C^J (u) is (E, m)-S } 



|G c( «>, u ,(u,v;£)| > e^™^^ 



< P 3£ e 



7 . 3 L o T c LQ (u J )(u j ,v j ;uj) ^ Q _ l{mMuN)Lo 



diBt(E )£ r(HW )) 

t i ' 



< P |3E G 7 : C^(u) is £~r} 

+ p{T Cio(Uj) K,, j;w ) > a-^e^"^^-^} 

T c l0 (« i )(wj,«j;w) 



<v 



3" w L" Ar e -'vC m ' i ''' JV ) I "'-' L o 



1/2 



< L- q + e -3°Lo 3 N L N ff/ ( m ,L ,N)L +L 1 ' 



-2p 



< 



l: 



-2p 



where we used property (WS1.A) (Wegner type bounds) to estimate the probability 
of i?-resonance in lines 4-5 and the bound (16.111) in line 6. □ 
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